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Abstract 

New examples of matrix quasi exactly solvable Schrodinger operators are 
constructed. One of them constitutes a matrix generalization of the quasi 
exactly solvable anharmonic oscillator, the corresponding invariant vector space 
is constructed explicitely. Also investigated are matrix generalizations of the 
Lame equation. 



1 Introduction 

The topic of Quasi Exactly Solvable (QES) operators [I], 0] has been developped in the 
last years. It consists of differential operators (mainly Schrodinger ones) which possess 
a finite dimensional invariant vector space, say V, of functions. So, the restriction of 
the eigenvalue equation to the space V leads to an algebraic eigenvalue problem. 



Although scalar QES operators have been classified in one || and several [|J 
variables, a classification of matrix QES operators is still missing. Yet, interesting 
examples of them have been obtained in relation with the stability analysis of soliton 
solutions occuring in some field theories |5], ^j. 

This problem was first addressed in and further developped in || and || . More 



recently H, 10], interesting tools for the classification of NxN matrix QES operators 



in one spatial dimension have been constructed and applied to the case N=2 |IT 
(although potentials of the Lame-type are not included in this classification). 

Here we consider a suitable class of finite dimensional vector spaces of N polyno- 
mials in a real variable and we construct families of operators preserving sub-classes 
of these vector spaces. The corresponding QES equations respectively constitute 
"coupled channel" generalizations of the scalar QES equations. 

Following the basic idea of QES operators |1J we consider the finite dimensional 
vector space of N-uples of polynomials of given degree rii, n 2 , ■ ■ ■ , in a real variable 
x. We slightly generalize this vector space by setting 

V = P 0P(m) © V(n 2 ) © ... © V{n N )) (1) 

where V(rii), % — 1, ...N denotes the set of real polynomials of degree at most in x 
while P is a fixed invertible NxN matrix operator; P can be interpreted as a change 
of basis in the vector space V[n\) © ... © V(n^). With such an interpretation, and 
if we assume n\ > n 2 > ■ • • > n^, it is reasonable to choose the matrix P as a lower 
triangular matrix with Pa = 1. 

2 QES Anharmonic matrix potentials 



We consider an operator of the form 



H(y) = --^K N + M 3 (y 2 ) (2) 
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where M 3 (x) is a iV x N hermitian matrix whose elements are polynomials of degree 
at most three in the argument x. After a standard "gauge transformation" of H(y) 
with the factor 

4>{y) = exp -{y2/ 4 + PiV 2 } (3) 

(pi,p 2 are arbitrary real parameters, p 2 > 0) and the change of variable x = y 2 , the 
operator equivalent to (|2|) reads 

, d 2 d s 

H ^ = -^ 2+2 Tx^ 

d 

+ 8{p 2 x + p^X—U-N 



(4pjx 3 + 8pip 2 x 2 + (4pi - 6p 2 )a; - 2pi)]I 



A' 



+ M 3 (x) (4) 

We now determine the form of the matrix M 3 (x) such that the operator H(x) possesses 
a finite dimensional invariant vector space of the type (P for generic values of pi,p 2 - 
In this purpose, the differential operators in the second line of (^) (i.e. 8p 2 x 2 d/dx and 
8p\xd/dx) have separately to be completed into operators which have this property. 
In order to make use of the results of || [Kj we conveniently rewrite H according to 

H{x) = -(4x-^ + 2-^)lhv 

dx z dx 

+ 8p 2 Q + + 8 Pl Q + 8p 2 {W 3 {x) + B) (5) 
where the operators Q + , Qo are defined by 

Q+ = x 2 ^- + 2xA-B , Q = x^- + A (6) 
and the constant matrices A, B are chosen to obey [A, B] = B, in such a way that 

[Q + ,Q ]=Q+- (7) 
Without loosing generality | JLU[| , we can choose A and -B in the form 

A = Diag(0,l,2,...,iV-l)-|lI Ar , B a , b = c b 5 a<b+1 (8) 



3 



The matrix W$(x) is defined by identification of (|j) and (|5]); it is symmetric and 
does not contain derivative. In the following, we further assume this matrix to be 
irreducible by ^-independent changes of basis; in particular, we exclude the cases 
where W 3 is diagonal. 



In H, [L0| , it was demonstrated that Q + ,Qo admit a finite-dimensional vector 
space; we now construct it explicitely. For this purpose, we define a family of vector 
spaces characterized by two integers N and p 

V(N,p) = V(p) ®V(p - 2) © . . . ®V{p - 2N + 2) (9) 

and an (invertible) N x N matrix operator P with matrix elements 

Pij = if i < j 
= 1 if i= j 

The vector space V = PV can be seen as a change of basis on the space V (for brevity, 
we do not write the dependence on p, N anylonger). 

The following proposition, which can be checked after an algebra, provides the 
explicit form of the invariant vector space of Q + , Q : 

Proposition 1 



p-^BP = P-\Q + + B)P-(Q + + B) (11) 

P- 1 Q + P = x 2 — + 2xA (12) 

ax 

P-'QoP = Qo (13) 

Because the operators on the right hand side of these equalities preserve V, it follows 
immediately that Q + and Q preserve the space V. 

The requirement that the (W3 + _B)-part of H also preserves V' is guaranteed by 
the Proposition 2 below. For later use we will note by J + , Jq, J_ the usual irreducible 
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representation of sl(2) by N x N matrices, i.e. 

[J , J±] = ±J± , [J+, J_] = 2J . (14) 
In particular we can set J = A — ^TrA 

Proposition 2 

If W%(x) is symmetric and irreducible, then the condition 

P~ 1 (W 3 + B)PVCV (15) 

is fulfilled if and only if 

B = cJ + , W 3 = -c(J+ + J_) (16) 

The values of the different parameters entering in the matrix B are therefore fixed 
up to the multiplicative factor c. 

So far, we proved that the condition is necessary and sufficient for iV up to seven 
and we are confident that it is true for arbitrary N. 

The final step is to check whether the operator in the first line of (|5|) also preserves 
V'. In this respect, it is sufficient to study the condition A-^V' C V. From the 
structure of the vector space V, it is quickly seen that this holds only for N=2. This 
shows that the generalization of the QES anharmonic oscillator to matrices is possible 
only for two-dimensional matrices and the 3-parameter family of operators of |J is 
recovered. 

Our results thus lead to a family of 3 parameters matrix quasi exactly solvable 
anharmonic oscillators (pi,p2,c). The parameters pi,P2 determine the anharmonic 
frequencies of the y 6 - and ?/ 4 -parts of the potential M 3 , while the constant c ensures 
a non-trivial coupling between the different oscillators. 
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3 Clebsch-Gordan coefficients 



The results presented in the previous section contains a formulation of the Clebsch- 
Gordan matrix of the tensor product of two representations of sl(2) in terms of dif- 
ferential operators. In this section we would like to further analyze some aspects of 
this realization. 

We consider on the one hand the standard representation of the algebra sl(2) 
by iV x iV-matrix generators noted J + , Jq, J_ and on the other hand the realization 
expressed in terms of differential operators : 

2 d d d 

j + = x — + 2/zx , jo = x— +fM , j- — — (17) 
dx dx dx 

With the choice \x = —n/2 (n integer), this realization preserves the vector space 
V(n) and constitutes the corner stone of QES equations 0. 

Then, we consider the tensor product of these representations, the generators of 
which are 

Q±,o = J±,o + J±,o • (18) 

Again, if /i = —n/2 , it is finite dimensional and preserves, in a reducible way, the 
vector space 

V = V(n) © V{n) ... © V(n) , N times (19) 

The explicit decomposition of the representation (|I5D, QI$D into irreducible ones is 
achieved in two steps. We define a first transformation by means of 

Q±,o = e xJ -Q±, e- xJ - , (20) 

leading to 

Q+ = x 2 -^- + 2x(J + utt) + J + , Q = x^- + J + /iH , Q- = (21) 

which can be identified with the operators (^) if A — Jq + /ill, i.e. if /x = N ~^~ 1 . 
In other words the two irreducible representations J e , j e under investigation have 
dimensions and p — N + 2, respectively. 



The second step leading to the desired decomposition reads 

Q±,o = P^QifiP (22) 

where the operator P is defined in eq.(|H]). The form of the operators Q is available 
from Proposition 1 : 

Q + =x 2 -^ + 2xA , Q = x-^ + A , Q- = ^ (23) 

It clearly reveals that the vector space (|9]) is preserved. As a consequence the matrix 
P cg = exp(— xJJ)P constitutes the Clebsch-Gordan matrix of the decomposition. 

It should be noticed that the operators J_, j_ preserve separately the vector space 
(|19D. Then, since the transformation (|20|) commutes with J_, it results that both J_ 
, j- also preserve V. This statement in fact provides the proof that the condition 
fll6|) of the Proposition 2 above is indeed sufficient. 



4 Lame- type operators. 

In the following the Jacobi elliptic functions 

sn(z, fe) , cn(^, k) , dn(z, k) . (24) 



of argument z and modulus k JT2[ are abreviated respectively by sn, cn, dn. These 
functions are periodic with period 4K(k),4K(k),2K(k) respectively {K(k) is the 
complete elliptic integral of the first type). 

It is well known that in order to reveal the algebraic properties of the Lame 
equation 

d 2 

+ k 2 N(N + l)sn 2 ]^(z) = Ei/)(z) , (25) 



L dz 2 

the relevant change of variable consists in posing x = sn 2 (z,k). In particular the 
second derivative operator is transformed into 

= Axil - x)(l - k 2 x)^- + 2(3k 2 x 2 - 2(1 + k 2 )x + l)^- (26) 
dz 2 d 2 x dx 
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and Eq. becomes a Fushs equation with four regular-singular points (at x = 

0,l/£; 2 ,l,oo). 

To our knowledge, attempts to construct (and classify) Schrodinger QES matrix 
operators with this type of change of variable have not been attempted so far. Par- 
ticular cases are emphasized in || The natural choice is to consider N x N 
Schrodinger matrix operators with potentials depending on the Jacobi ellipitc func- 
tions and which possess algebraic properties similar to Eq. (|25|) . More specifically, we 
consider operators of the form 

H(z) = + V D (z) + V T (z) (27) 

with 

V D = sn 2 diag(ai, a 2 , • • • , a N ) + diag(&i, b 2 , . . . , b N ) (28) 

where aj, bj denote real constants (without loosing generality we assume J2f=i bj = 0) 
and Vj is a symmetric off-diagonal matrix of the form 



9ij sn Ql ^cn a2 ^'dn a3lJ if i ^ j 
if % = j 

Owing the periodicity of the Jacobi elliptic functions, the family hamiltonian above 
is to be considered on the Hilbert space of periodic functions on [0,AK(k)]. 

Unfortunately, it has not been possible to classify the QES operators of the form 
(|27|), here we will describe the cases that we were able to treat. The difficulty of 
achieving the classification will appear from these few examples. The properties of 
the Jacobi functions that are useful to make the calculations are listed in the appendix. 

4.1 case N=l 

This case corresponds to the Lame equation which was discussed lengthly e.g. in |f[4}| . 
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4.2 Case N=2 

The non-diagonal potential reduces to only one component 



Vi = 6 sn ai cn Q2 dn Q3 (7i (29) 

Using a similarity transformation of the form 

H(x) = U- 1 (z)H(z)U(z) , U(z) = diag(sn /31 cn' 32 dn & ,sn 71 cn 72 dn 73 ) (30) 

sets the operator ( |2"TD into a form with polynomial coefficients in the variable x = sn 2 
provided (see Appendix) 

• Phli = or 1 , j = 1,2,3 

• atj ± (Pj — 7j) = non-negative even integer , j = 1, 2, 3. 

Considering first the cases where V\ 2 is a constant (01,2,3 = 0) or a single linear factor 
in one of the Jacobi functions, we were able to show that the only possible QES 
operators are available for 9 = 0. 

For the three cases a) Vu = 9 sn cn, b) V12 = sn dn and c) Vu = 9 cn dn, the 
construction of non-decoupled QES operators is possible. One remarkable feature is 
that in each case, two sets of values of the coupling constants ai,a2,b±,9 lead to four 
different algebraizations of the corresponding operator. We now discuss it in detail 
(posing b = b\ = —b 2 ). 

4.2.1 Case V 12 = 9 sn cn 

The two sets of values of the coupling constants leading to QES operators are noted 
Type 1 and Type 2. 

Type 1 

Oj = k 2 (4m 2 + 2m + 1) - 2b 
a 2 = k 2 (4m 2 + 2m + 1) + 2b 
9 2 = Ab 2 - k\l + Am) 2 
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The parameters b and k remain free, m is an integer. Four invariant spaces are 
available : 



Vi 



V 2 



V, 



V 4 = 




sn cn dn 
dn 

Type 2 

ai = k 2 (Am 2 + 6m + 1) - 2b 
a 2 = k 2 (Am 2 + Qm + 1) + 2b 
e 2 = Ab 2 - k 4 {3 + Am) 2 

The corresponding invariant vector spaces read 



-9 



k x = 



26 + k 2 (1 + 4m) 
9 

2b- k 2 (l + 4m) 
1 



K-3 = ~- 



K-2 



K4 = 



(31) 

(32) 
(33) 
(34) 




V, = 



v fi 



V 7 = 



V 8 = 



4.2.2 Case V 12 = 9 sn dn 



k 5 = 



k 2 {Am + 2,) -2b 
9 



k 2 (Am + 3) -2b 



,K 7 = 



9 

2b-k 2 (Am + 3 
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(35) 
(36) 
37) 



In the following cases we just write the values of the coupling constants leading to 
QES operators; the corresponding invariant vector spaces are similar to the ones 
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above, their explicit form can be obtained in a straightforward way. 
Type 1 

Ql = k 2 (4m 2 + 2m + 1 - 2b) 
a 2 = k 2 (Am 2 + 2m + 1 + 2b) 
6 2 = 4k 2 b 2 -k 2 (l + 4m) 2 

Type 2 

ai = k 2 (4m 2 + Qm + 1 - 2b) 
a 2 = k 2 (4m 2 + Qm + 1 + 2b) 
9 2 = 4k 2 b 2 -k 2 (3 + 4m) 2 



4.2.3 Case V 12 = 6 cn dn 
Type 1 

ai = k 2 (4m 2 + 2m + 1) - 2bj^ 
a 2 = k 2 (Am 2 + 2m + 1) + 26^ 
# 2 = A; 2 (l + 4m) 2 -46 2 T |^ 

Type 2 

ai = A; 2 (4m 2 + Qm + 1) - 26^ 
a 2 = k 2 (4m 2 + 6m + 1) + 26^ 
# 2 = P(l + 4m) 2 -4& 2 T ^ 

Concerning the operators, we want to point out the following things: 

• The extraction of the prefactor is done in two steps. After factorizing the 
appropriate products of Jacobi functions, the non-diagonal part of the potential 
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takes the form 

x 2 ~ e 



2 or 1 or (39) 



x e 

which is clearly incompatible with an operator preserving a vector space of 
the form V(n) © V(m) for integers m, n. The setting of the operator in a form 
preserving such a vector space is realized by the second prefactor (the triangular 
matrix one). As an example if e = in (|39|), then use is made of the relations 



(40) 



(41) 

with A x in order to eliminate the quadratic power of x occuring in V\2- 

For V\2 = sn cn, the limit k — > is non trivial and leads to a potential of the 
form 




cos 2 (x) — I cos(x)sin(x) 



(42) 



\ cos(x)sin(x) sin 2 (x) — \ 
In particular, the dependence on m disappears, so this integer can be choosen 
arbitrarily large and an infinite set of algebraic eigenvectors occurs in this limit. 
This confirms the fact that the corresponding Schrodinger equation is com- 
pletely solvable and is related to the stability of some soliton solutions occuring 
in the Goldstone model in 1+1 dimensions with a periodic condition for the 
space coordinate. 

The case V%2 = 9 cn dn with b = (1 + k 2 )/2 was studied lengthly in p. The 
corresponding Schrodinger equation is related to the stability analysis of the 
sphaleron solution in the U(l)-Abelian Higgs model in 1+1 dimensions (again 
with a periodic condition for the space coordinate). 
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4.3 Case N=3 



We attempted to construct 3x3 matrix operators with some choice for the constants 



a 



aij 



(a,i,j = 1,2,3) inspired from the results above. Several trials were unsuccesfull 



(namely with V i2 -r- V 23 -r- sn cn); however the choice 

/ 



V 



Q\ cn dn V23 
9i cn dn 9 2 cn dn 
V 23 9 2 cn dn 



(43) 



leads to the wanted form of equations and preserves a vector space of the form 

/ __ n n \ 



V = 



The condition 



cn 











dn 











cn 



1 1 a 
1 7 
1 



if(z)V CV 



( V{n-2) X 
7>(n- 1) 
V{n) 



(44) 



(45) 



leads to ten equations for the parameters a, (3, 7, fc, aj, b; L and These equations are 
compatible with each other only if k — 1 and V23 = ^cn 2 . We found it convenient 
to leave the parameters a,/3, 7 free and to express the coupling constants cii,bi and 
0a in terms of them. We refrain to write these tedious expressions in general but we 
mention that they have the polynomial 



(1 + 7 2 + /9 2 )(1 + P 2 + 7 2 + «V - 2^7) 
as common denominator. Also they are such that 

d\ + a 2 + as = 2(6n 2 — 3n + 4) 



(46) 



(47) 



For the case with e.g. a = = 7 = 1 we find 

12n 2 - 10n + ll 



ai= a 2 = 



, a 3 = -(6n 2 + n + 1) 



4n-3 

bi = b 2 = — - — , 3 = -2bi 



(48) 
(49) 
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e 1 = l^ , 92 = -iH±i , fe ^ 2(4 " + 1) (50) 

6 ' 3 ' 3 v ; 

as another example, if a = —f3 = 1 , 7 = 0, we have 

2(6n 2 -7n + 3) 1/9 
ai = - 3 , a 2 = a 3 = -(12n 2 -2n + 9) (51) 

^ = ^1) , , * = (52) 

* = ^ , fc = 2^ , *3 = ^^ (53) 
3d 3 

Remarkably, replacing the first factor in (fP|) by diag(dn, cn, dn), we obtain the 
same solution for ai,a 2 , .... So the operator ( f43|) possesses at least a double alge- 
braization, but we have not attempted to construct the other possibly existing ones 
yet. 



The 3x3-matrix QES potential (43) strongly contrasts with the 2x2 ones obtained 
above. While the former depends on three free parameters but exists only on the full 
line (since k — 1), the latter can have an arbitrary period (equal to 4K(k)) but 
has only one free parameter (for instance noted by b) for a fixed k. The fact that 
the operator H(z) for a 3x3-matrix QES potential is quasi exactly solvable only for 
k 2 = 1 was unexpected to us. Accordingly, this operator corresponds to a 3 x 3 matrix 
version of the Poschl- Teller operator. 



5 Concluding remarks 

In the first part of this paper we obtained the matrix generalization of the celebrated 
sextic QES anharmonic oscillator. The examples of operators presented in the second 
part reflect the difficulty to classify the coupled-channel (or matrix) QES Schrodinger 
equations when the change of variable involves several singular points like in (|2Tj|). 
We hope that this note will motivate further investigation of this problem. 
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6 Appendix A 



For k = we have K(0) = | and the Jacobi functions reduce to standard trigono- 
metric functions : sa(z, 0) = sin(z), cn(z, 0) = cos(z), dn(z, 0) = 1. For = 1, we 
have if (1) = oo and the Jacobi functions reduce to elementary functions : 

sn(z, 1) = tanh(z) (54) 

™{z, 1) = dn(z, 1) = ^— (55) 

In this limit, the Lame equation becomes a Poschl- Teller equation. For generic values 
of k, the Jacobi functions obey the following relations : 

cn 2 + sn 2 = 1 , dn 2 + A; 2 sn 2 = 1 (56) 

— sn = cn dn , — cn = — sn dn , — dn = — k 2 sn cn (57) 

dz dz dz 

These identities as well as the following ones are useful to establish the equations 
in the variable x = sn 2 after the prefactor including the Jacobi functions has been 
extracted : 

/ /"// (sn cn dn) /'// 

1 

sn 2k 2 x - (1 + k 2 ) k 2 x 2 - (1 + k 2 )x + 1 
cn 2k 2 x - 1 k 2 x 2 - x 

dn 2k 2 x - k 2 k 2 x 2 - k 2 x 
cn dn 6A; 2 a; - (1 + k 2 ) 2k 2 x 2 - (1 + k 2 )x 
sndn Qk 2 x- (l+Ak 2 ) 2k 2 x 2 - (1 + 2k 2 )x + 1 
sn cn Qk 2 x - (4 + k 2 ) 2k 2 x 2 - (2 + k 2 )x + 1 
sncndn 12k 2 x - 4(1 + k 2 ) 3k 2 x 2 - 2(1 + k 2 )x + 1 

No similar identities are available (to our knowledge) with different choices of the 
function /. 
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